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Abstract 

We show that the recent string theory embedding of a spacetime with nonrelativistic 
Schrodinger symmetry can be generahsed to a twenty one dimensional family of solutions 
with that symmetry. Our solutions include IIB backgrounds with no three form flux 
turned on, and arise as near horizon limits of branewave spacetimes. We show that 
there is a hypersurface in the space of these theories where an instability appears in 
the gravitational description, indicating a phase transition in the nonrelativistic field 
theory dual. We also present simple embeddings of duals for nonrelativistic critical 
points where the dynamical critical exponent can take many values z ^ 2. 



1 Background 



In recent studies of the (strongly coupled) quark-gluon plasma, just above the QCD de- 
confinement temperature, the AdS/CFT correspondence [1] is finally starting to deliver 
qualitative and even quantitative results of experimental relevance that are difficult to ob- 
tain by more traditional theoretical methods (for overviews see for instance [2, 3, 4]). 

The best understood versions of the AdS/CFT correspondence describe a duality be- 
tween a conformally invariant gauge theory and string theory in an asymptotically Anti-de 
Sitter space background. Taking a large N and (if necessary) strong coupling limit in the 
gauge theory corresponds to taking the Anti-de Sitter space to be weakly curved, and hence 
string theory may be described by a ten or eleven dimensional supcrgravity theory. 

QCD is a gauge theory, and so it is natural to hope that certain quantities will be 
in the same 'universality class' as the supcrsymmetric gauge theories that can be studied 
using AdS/CFT. This hope seems to be bourn out for the high temperature deconfined 
phase of QCD [2, 3, 4]. On the other hand, there are many quantities of interest that are 
sensitive to the fact that QCD is not a conformal theory, but rather asymptotically free and 
confining. Extensions of the AdS/CFT correspondence to non-conformal theories are well 
studied. However, they generically suffer from the feature that parametrically decoupling 
the confining physics from other scales in the string background (such as the Kaluza-Klein 
scale) requires going beyond the supcrgravity approximation. This problem is reviewed for 
instance in [5]. 

One is then led to wonder whether there are exactly conformally (or at least, scale) 
invariant systems in nature that would be kinematically closer to the simplest versions of 
AdS/CFT. Conformal invariance has in fact been observed to arise in an increasing number 
of condensed matter systems - in the vicinity of a quantum critical point [6, 7, 8]. 

The microscopic description of these condensed matter systems is certainly not relativis- 
tic. It is remarkable therefore that certain quantum critical points exhibit an emergence 
of the full relativistic conformal symmetry. This occurs for instance in certain superfluid- 
insulator transitions and at certain 'deconfined' quantum critical points (see for instance 
[9, 10] for references). If a nearby relativistic quantum critical point is believed to play a 
role in some dynamical process, then one can hope that standard AdS/CFT techniques can 
be applied. This was the methodological assumption in [10, 11, 12], which studied certain 
transport processes in unconventional superconductors. 

It is also common for quantum critical points to display a non-relativistic scale invari- 
ance. These theories have a scale invariance that treats space and time anisotropically, as 
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captured by the dynamical critical exponent z. Thus for instance, as the critical point is 

approached, the mass gap A, which determines relaxation timescales, and the coherence 
length ^ are related by 



(1) 

The case z = \ corresponds to a relativistic theory. For general z the symmetry of the 
theory is simply the Galilean group enhanced by a dilatation operation. The case z = 2 
is also singled out, as the algebra may be enhanced in this case by an additional 'special 
conformal' generator. This is called the Schrodinger algebra and enforces strong constraints 
on the dynamics. We will not review this algebra in detail, for recent discussions see for 
instance [13, 14, 15]. Early discussion of the symmetry of the Schrodinger equation are 
[16, 17] and an early geometric discussion of these symmetries appears in [18]. Examples of 
experimental systems with z = 2 include fermions at unitarity, see [14, 15] for references. 

Recent work has suggested that the AdS / GET correspondence can be adapted to study 
strongly coupled non-relativistic quantum critical points. We will concentrate in this paper 
on the most constrained case of z = 2, and on field theory in 2 + 1 dimensions. We shall 
also present some interesting extensions of our results to general z. The starting point is to 
write down a metric upon which the Schrodinger symmetry is geometrically realised. The 
following five dimensional metric does the job [14, 15] 



(2) 



Here o" is a parameter measuring deformation of the metric away from Anti-de Sitter space. 
Furthermore, one should periodically identify 

x~ ^ x~ + 27rr~ . (3) 

This allows P_ = idr^.- to be interpreted as particle number 

with an integer. If wc then consider .t+ to be time, so that P+ = 18^+ is energy, then 
positivity of energy for classical fields in the spacctime (2) requires A > 0. 

The works [14, 15] showed that by generalising standard AdS/CFT arguments to fields 
propagating in the background metric (2) one obtains sensible results for correlators of op- 
erators in a Schrodinger invariant theory. In order to exploit the full power of the AdS/CFT 
correspondence, and obtain, for instance, the finite temperature and number density back- 
grounds, it is necessary to embed the metric (2) into a consistent theory of gravity, such as 
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string or M theory. This was achieved in the works [19, 20, 21] who obtained the following 
solution of IIB string theory on AdS^ x X^, with a Sasaki-Einstein space 

dsjiB = (^r^ {-2dx+dx- + dx^) + ^ - a\\dx+f + ds\_E + > (5) 
B2 = ae^/^R^r^dx+ Ar), (6) 
F5 = 4i?4(l+*)vol(X5), (7) 

where $ is the dilaton, which is constant. In this expression the Sasaki-Einstein space has 
been written as a fibration over a local Kahler-Einstein manifold 

ds\, = dsl_E + ri" , (8) 

such that df] = 2lo, the Kahlcr form on the Kahler-Einstein base. In the maximally sym- 
metric case X^ = S^, the Kahler-Einstein manifold is CP^. Closely related solutions had 
been found previously in [22, 23]. 

Before moving on, two comments are in order. Firstly, it was emphasised in [20] that 
the circle that wc have periodically identified in (3) is null and therefore has zero proper 
length. This means that strings winding the circle are light and strictly speaking one cannot 
trust the supergravity regime. It was also noted how this situation can be ameliorated by 
considering finite density, as is indeed appropriate for comparison to experimental setups. 
Secondly, an alternative proposal for realising a gravity dual of non-relativistic theories was 
presented in [24, 25]. These papers use pure AdS, with the relativistic symmetry being 
broken by boundary conditions. This approach has been critiqued, in our view correctly, in 
[20, 19]. 

In this work we shall present some new solutions to IIB supergravity with a Schrodinger 
symmetry. There will be a large family of solutions (21 dimensional in the case of AdS^ x S^) 
that interpolate from the solution (5) to a solution which only involves the metric. As we 
move along the space of solutions, an instability develops on a critical hypersurface in 
parameter space. This instability is presumably dual to a (quantum) phase transition in 
the coupling space of the quantum field theories. 

The single step in our construction of these solutions is to start with AdS^ x X^ and then 
to introduce a harmonic function on X5 into the AdS^ part of the metric. An intriguing 
feature of this construction is that the resulting exponent z is directly given by an eigenvalue 
of the Laplacian on X5. In general z ^ 2. We go on to show how such spaces arise as the 
near horizon limit of a branewave traveling along N D3 branes. We also show that despite 
the mixing of X5 and AdS^ coordinates, a scalar field propagating in the ten dimensional 
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background is equivalent (in the z = 2 case) to a scalar field in the five dimensional space 

(2) but with a renormalised mass. 

Towards the end of this paper we construct two types of explicit renormalisation group 
flow solutions that are flows away from an ultraviolet (UV) Schrodinger invariancc. The first 
is a novel solution that appears to describe our field theories with a finite energy density. 
The second are nonrelativistic deformations of the known Coulomb branch solutions of 
A/" = 4 super Yang-Mills theory. We conclude with comments on the dual field theory 
interpretation of our solutions, and a rather long to do list. 



2 A family of IIB Schrodinger solutions 
2.1 Solutions with no fluxes 

Consider the following IIB ansatz, which has the full Schrodinger symmetry 

ds]iB = [r^ {-2dx+dx- + dx^) + ^ - f{X5)r\dx+f + dsj^^^ , (9) 
F5 = 4i?^(l+*)vol(X5). (10) 

Here X5 can be any five dimensional Einstein manifold, i.e. = ^g^- What we have 

done is to insert a function /(X5), depending on the coordinates on X5, into the noncompact 
part of the metric. Thus the spacetime describes a wave propagating in AdS^ with a profile 
along the internal directions. A similar ansatz was considered in [26]. In fact, section 2.2 
of their paper describes a specific example of the solutions we shall discuss in the following 
subsection. 

The wave in AdS^ in (9) only alters the R++ component of the Ricci tensor relative to 
the pure AdS^ background. It is easy to check that our ansatz solves the Einstein frame 
IIB equations of motion 

Rmn = —FmabcdFn^^''^^ , (11) 

provided that 

-^xj = m- (12) 

If / diverges at any point this will introduce singularities into the spacetime. These are not 
curvature singularities, but rather infinite tidal forces that appear for geodesic observers. 
For a recent discussion see for instance [27]. This is sufficient for time evolution in the 
spacetime to be ill-defined, so we should require that / is regular. Therefore (12) implies 
that / should be a scalar harmonic on X5 with a specific eigenvalue. 
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When we consider general z below, we will find that the scalar harmonic eigenvalue 
depends on z. For the moment we can ask which Einstein manifolds have a scalar harmonic 
with eigenvalue 12, as required by (12). For we have the spectrum 

-vis = lit + 4) . (13) 

So the 1 = 2 harmonics give solutions. The degeneracy of £ = 2 on 5^ is 20, so AdS^ x has 
a twenty dimensional family of solutions of the form (9). The degeneracy can be obtained 
as the number of tracefree quadratic forms in six variables. 

In fact, any Sasaki-Einstein manifold has eigenvalues of the form (13), see for instance 
[28]. However, in general t does not run over all the positive integers, so we are not 
guaranteed the existence of a mode with eigenvalue 12. A familiar case where the spectrum 
is again known exactly is the homogeneous Sasaki-Einstein manifold T^'^. The spectrum is 
(see e.g. [29]) 

-V|m = 6 + 1) + h{h + 1) - . (14) 

Here r is an integer, and £1,^2 can be integers or half integers, constrained as specified in 
for instance [29, 30]. There are precisely two combinations with the required eigenvalue: 
(r, £1,^2) = (0,1,0) and (0,0,1). Because i\ and I2 label standard spherical harmonics on 
S^, this means that we have a total of six modes in this case, and hence a six parameter 
family of solutions. 

It is interesting to note that if an Einstein manifold admits eigenvalues with — = 12, 
then the corresponding AdSf, x X5 compactification has a scalar mode in AdS^ which 
saturates the Breitenlohner-Preedman bound for stability. This is because certain metric 
and five form fluctuations about AdS^ x X5 have masses m^i?^ = — + 16ib8y^— V^^ + 4 
(see e.g. [31]). Putting — = 12 leads to the Breitenlohner-Freedman bound for AdS^: 
n?^ = -4. 

On any X5, a nontrivial spherical harmonic has regions in which it is both positive 
and negative. This follows from orthogonality of the harmonic with a constant function. 
Therefore our function f{X^) appearing in the spacetime (9) will be negative in certain 
regions of the internal space. One might therefore worry about stability [27, 33] and geodesic 
completeness [32] of the spacetime. We will return to these points below. First we shall 
show how NSNS flux can be added to our solution (9) to allow an interpolation between 
these solutions without flux and the known solutions (5). 
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2.2 Solutions with fluxes 

We can now add a NSNS two form potential to our solution (9) . This will allow the family 
of solutions to connect onto the solutions (5). The ansatz is now 

dsjiB = [r^ {-2dx+dx- + dx^) + ^ - f{Xr,y{dx+f + ds]i_E + rf^ , (15) 
B2 = (7e*/2i?Vdar+ A r? , (16) 
F5 = 4i?4(l+*)vol(X5). (17) 

Here we again require X5 to be Sasaki-Einstein and use the same notations as above. 

It is easy to check that the five form and three form field equations are satisfied exactly 
as they were for the solution (5), the function /{X^) does not appear in these equations. 
The Einstein equation becomes 

Rmn = -^FmabcdFn^^^^ H — (^mabHn^^ — J^gMNHABcH'^^'^^ , (18) 

where = dB2 and $ is the (constant) dilaton. The last term in this equation vanishes 
because is null in (15). The Einstein equation is solved provided that 



Prom this equation we obtain 



-VL/ = 12/-12(7^ (19) 



f = cT^ + f, (20) 



where / is, as previously in (12), a scalar harmonic on X5 with eigenvalue 12. We see that 
these solutions are precisely a linear superposition of our solution (9) without fluxes, and 
the known solution (5). By tuning a and the magnitude of / we interpolate between these 
two limits. Thus we obtain a 21 dimensional family of solutions from AdS^ x and a 7 
dimensional family from AdS^ x T^'^. We can also note that by taking a to be sufficiently 
large relative to /, the function / appearing in the metric will become positive everywhere. 
This fact will shortly have consequences in our stability analysis. 



3 Scalar field fiuctuations 

One might have expected that the mixing of internal and noncompact coordinates in our 
solutions (15) would complicate the study of fluctuations about the background. This is 
not the case. We will now show that, for scalar fields at least, the only effect of the mixing 
is to 'renormalise' the effective Kaluza-Klein mass of the modes. 
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The linearised equation of motion for a (massive) ten dimensional scalar field is 

V^(/) = m^(^. (21) 
It is easily seen that we can solve equation (21) by writing 

(j) = ^r)Yx{X5) e-''^^++^'=-="-^A^x- _ (22) 

Here oj is the frequency, k the wavevector and M the rest mass. The eigenvalue A is the 
solution to a Schrodinger-like equation on 

[-V\^+M^f]Yx = \Yx. (23) 

The function $(r) solves 

cf^ 5 2Mu - fc^ ^ _ + 

dr'^ r dr 

Therefore we see that the only effect of the mixing of internal and noncompact coordinates 
is to 'renormalise' the effective mass of the five dimensional field: — > + A. The wave 
equation separates even though the background is not a direct product. If the rest mass 
M = 0, then there is no dependence on / and Y\ arc simply the spherical harmonics on X^. 
For nonzero M, the functions Yx solve a modified equation on X^. 

The five dimensional equation (24) is easily solved [14, 15] to give, up to normalisation, 
the modified Bessel function 

,25, 

where 

p=\fW^^2M^, V = + A + 4 . (26) 

We have imposed regularity of the solution at the horizon r = (assuming that p is real). 
The different modes of this equation correspond to scalar operators O with scaling dimension 

Ao = 2 + i.. (27) 

If < < 1 we can also take = 2 — v [14]. Thus we see explicitly that the effect of our 
function / is to shift the scaling dimension of the operators dual to bulk scalar fields. 

4 Stability and geodesic completeness: a phase transition 

We noted that in our solution without fiux (9), the function / appearing in the metric would 
necessarily be negative over some region of X^. This is potentially associated with various 
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spacetime pathologies, see e.g. [27, 33], two of which we shall now investigate. On the other 
hand, in the previously known background (5), / is positive everywhere. In our family 
of solutions (15) interpolating between these limits, / becomes negative on some critical 
hypersurface in the space of solutions. A pathology arising on this surface would indicate 
an interesting zero temperature phase transition in the dual family of nonrclativistic field 
theories. We will see that while the space remains geodesically complete across the transition 
to nonpositive /, there is indeed an instability once / becomes sufficiently negative. 

4.1 Appearance of an unstable scalar fluctuation 

The first question we consider is stability against scalar field perturbations, such as the 
dilaton. Stability in pp wave backgrounds is complicated by various factors such as the lack 
of global hyperbolicity [34]. See however the comments in, for instance, [33]. We shall take 
the following simple approach that seems appropriate for the present context. 

Firstly, we fix the sign of the rest mass M > 0. This is a choice of direction of the 
Killing vector generating translations in x~ . Next, perform a simple change of variables on 
the effective five dimensional equation for a scalar field in (24) above. The objective is to 
put the equation in Schrodinger form. Let 

u = -, $ = u^/^^ . (28) 
r 

Furthermore, we will put = 0, as zero momentum gives the potentially most unstable 
mode. Then from (24) we obtain 

Recall that v is given by (26) above. This is a very well studied Schrodinger equation (with 
'Energy' E = 2Muj). In particular, it is easy to see that when i/^ < there is a continuum of 
bound states with arbitrarily low energies. This means that there are normalisable solutions 
to (29) with arbitrarily negative Muj. 

An instability is usually signalled in general relativity by an exponentially growing mode, 
in which uj is imaginary. This is not what we have here. In the usual case of pure AdS 
space, / = 0, the solutions we have found for i/^ < correspond to modes with masses below 
the Breitenlohner-Freedman bound [35], that are growing exponentially in the boundary 
Minkowski time t = x'^ + . In the present setup, we have periodically identified x^, and 
hence the rest mass M must be real. This constraint, and self-adjointness of the Schrodinger 
equation, requires u to be real also and thus there is no exponential growth. 
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However, given that we have taken M > 0, then Mu negative requires oj negative. We 
will see shortly that large M is necessary to obtain h'^ < 0. Nonetheless, we can take M large 
but finite. The fact that Mlo is unbounded below then implies that we have normalisable 
solutions with u arbitrarily negative. The frequency cu is nothing but the energy of the dual 
field theory, with respect to the nonrelativistic time . Our analysis will therefore imply 
that the Hamiltonian of the dual nonrelativistic theory becomes unbounded below when 

< 0. This is indeed an instability. We do not see exponential growth because we have 
taken a free scalar field. Generic interactions will send the system into energetic free-fall. 

We can also note from (27) above that u"^ < corresponds to an imaginary scaling 
dimension for the dual operator O. This reinforces the idea that the system is unstable in 
these cases. 

The question is therefore: when do modes appear with i/^ < 0? Recall that = 
+ A + 4. The standard AdS^ x Xr, vacua of IIB string theory are stable provided that 
there are no five dimensional modes with masses below the Breitenlohner-Preedman bound: 
+ 4 > 0. Let us assume we are compactifying on such an X5. This class of five manifolds 
includes all Sasaki-Einstein manifolds. Our solutions will therefore become unstable if A 
can become sufficiently negative. Recall that the allowed values of A are eigenvalues of the 
operator (23) on X5. It is clear that if / is positive in (23) then A must also be positive. 
Thus 

/>0 => Stability (30) 

Of course we have only shown stability against the particular mode (21). 

We will now show that if we are sufficiently close in parameter space to the solution 
(9), with no three form fiux, then by taking M sufficiently large an instability will appear. 
We will not show this in generality, but rather for the case of S^. Let us scale M to be 
parametrically large, which we can do, as the theory contains excitations with arbitrarily 
large particle number. It then follows from (23) that a sufficient condition for an arbitrarily 
negative A is if we can find a test function Y on such that the potential energy is negative. 
That is 

3Y s.t. / fY^dn < Instability. (31) 

The potential term will eventually dominate over the kinetic term in (23) if we take M to 
be sufficiently large. 

If we introduce cartesian coordinates {x^} on M^, then we have from above that 

f = a^ + ttijx'x^ , (32) 
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with Qij a symmetric tracefree real matrix. The function / is of course restricted to the 

unit S'-' C M*^. Let us now take as a test function an i = 1 harmonic on the five sphere. 
That is: Y = h^x^ . This gives 

/ fY'^dn = a'^hkhi{x^x^)s5 + aijbkbi{x''x^x^x^)s5 

= V„l(S=)(^!Mi + ^). (33) 

Denote the eigenvalues of a-ij by a^'^\ A tracefree symmetric matrix necessarily has at least 
one negative eigenvalue. Therefore 

4(7^ < Imina^l ^ Instability. (34) 

i 

This relationship shows that if a (i.e. the three form flux) is sufficiently small, the theory 
becomes unstable, as advertised. It may well be that the theory becomes unstable at a 
larger value of a than that given in (34). The results (30) and (34) leave the stability in 



the range G 



^ I miuj a^*^ I , I miuj a^*^ | 



undetermined. 

To summarise the results from this section: Equations (30) and (34) give conditions 
on the function appearing in (32) to give stable and unstable spacetimes, respectively. We 
have not found the precise hypersurface in {cr^,ajj} space where an instability sets in, but 
we have bounded it on both sides. 

It is worth highlighting that this instability is catastrophic for the spacetime. As soon as 
the the criterion a < cicrit. (fs^*-*) is satisfied then all modes in the theory with sufficiently large 
particle number become unstable at once. It would be extremely interesting to understand 
the field theoretic dual description of this process. 

4.2 Geodesic completeness 

Another potential pathology of the spacetime that can arise when / is not positive is geodesic 
incompleteness. In particular, we have to check that no timelike or null geodesies can reach 
infinity in a finite affinely parametrised worldline time [32] . For simplicity let us consider a 
geodesic that is at a constant location on X5. This is possible at the stationary points of 
/(X5). It is easy to check that relaxing this restriction will not alter our conclusions - the 
geodesic would simply oscillate in the X5 directions as it moves out towards infinity. 

It is straightforward to check that to leading order at large r, (timelike or null) geodesies 
satisfy 

-{M'^ + fE'^y, (35) 



ds 
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where s is an afHne parameter and and E"^ are positive constants (M can be zero) . If / 
is positive everywhere, then this equation is inconsistent and hence geodesies cannot reach 
infinity at all. If / is negative then, by taking E sufficiently large, the right hand side of 
(35) can become positive. The equation is easily solved to give 

r ~ 6^"% (36) 

for some positive constant m. Therefore we see that it takes an infinite affine time for 
geodesies to reach infinity. The spacetime is geodesically complete in all cases. This result 
is closely related to the regularity of plane wave spacetimes [36] . We shall see the connection 
to plane waves explicitly in section 6 below. 



5 Many solutions with z ^ 2 

Generalising our above construction to a large number of values of z is straightforward and 
leads to an interesting observation. Consider the ansatz 

^iiB = [r'' {-2dx+dx- + dx^) + ^ - f{X^y'{dx+f + ds],_E + rf^ , (37) 
Fs = 4i?4(l+*)vol(X5). (38) 



Note that we have not included a B2 field. The B2 equation is not easy to satisfy for general 
z, which is one of the reasons why the papers [19, 20, 21] concentrated on the case z = 2. 
With the above ansatz (37) the Einstein equation is satisfied provided that 

-V^5/ = 4(^2-1)/. (39) 

Once again, we have regular solutions whenever the manifold admits a scalar harmonic 
with a specific eigenvalue. 

For the AdS^ x case, equation (39), together with (13), implies that we have solutions 
for all 

(40) 

That is to say, z can take an integer or half-integer value greater than one. The dimension- 
ality of the family of solutions obtained is given by the number of tracefree rank £ symmetric 
tensors on M^: 

m = ^{l+£){2+ ef{3 + £) = t^^f_zA . (41) 
Clearly this becomes very large for large z. 
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For a general X5, the allowed values of z will be a sequence of irrational numbers, di- 
rectly related to the scalar eigenvalues of X5 through (39). This gives a rather intriguing 
connection between spectral properties of Einstein manifolds and dynamical critical expo- 
nents. A general statement that can be made is the following. The Lichnerowicz theorem, 
see e.g. [37], applied to five dimensions, says that the lowest nontrivial Laplacian eigenvalue 
of a five dimensional Einstein manifold with i?^^ = ^g^n bounded below 

-Vl,,>5, (42) 

with equality only for the case of . Using our relation (39) we obtain 

z>\. (43) 

Of course there is also the z = \ possibility of pure AdS which corresponds to the triv- 
ial (zero) eigenvalue. Thus our construction will only enable us to obtain theories with 
dynamical critical exponents satisfying (43). 

We should note the following two facts. Firstly, for z > 2 the spacctimc (37) is not 
geodcsically complete. It is easy to sec that at large r, and in regions of X5 where / is 
negative, there are (null and timelike) geodesies that behave like 

" (S - So)V(--2) • 

Again, s is an affine parameter and sq a constant. These geodesies reach infinity at a finite 
affine parameter s = sq. The geodesic incompleteness may possibly be related to the lack 
of scaling near infinity for z > 2 that was noted in [15]. It is not clear at this stage how 
important this result is given that the large radius region is problematic in any case [20]. 
This is because, even with a putative solution at finite number density, at large radius the 
periodic identification of the x" circle means that wound strings become arbitrarily light 
and the supergravity description breaks down. In our solution of course, without number 
density, the circle is null everywhere. 

Secondly, if z < 2 the spacetime is singular at the origin r = 0. Although the curvatures 
are finite, there are pp singularities in regions where / is negative. For instance, a radially 
infalling null geodesic experiences infinite tidal forces as r ^ in a parallelly propagating 
orthonornial frame. These are mild null singularities that may well disappear in a finite 
temperature solution. 

We will not study these solutions further here. It would be interesting to investigate 
their stability. Our experience in the z = 2 case above and the absence of fiux in these 
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solutions might suggest that they are unstable. However, unlike we found in section 3 
above, the ten dimensional wave equation in the background (37) does not separate nicely 
for z ^ 2, which complicates the analysis. 



6 Undoing the near horizon limit: branewaves 

It is instructive to write down a IIB geometry that asymptotes to a pp wave in flat space 
and which develops a Schrodinger symmetry in a near horizon limit. The solution is 

ds]iB = {-^dx+dx- + da;2 - f{X5)r\dx+f) + H{rfl'' {dr^ + rV^,) (45) 

ds\i = ds\_E + rf, (46) 

B2 = ae*/V2dx+ A ry , (47) 
—H'(r) 

F5 = -rrr-y-il + *)dr A dx^ A dx~ A dx^ A dx^ . (48) 

Here H{r) is the usual extremal D3 brane function 

H{r) = 1 + ^ . (49) 

It is straightforward to check that the equation for the three form field strength is indepen- 
dent of the function H, and therefore remains satisfied for this more general metric ansatz. 
The Einstein equations (18) are satisfied provided that, once again, 

-Vj,J = 12f-12a\ (50) 

As above, this is solved by letting f = cr^ + f, with / being a scalar harmonic on with 
eigenvalue 12. The near horizon limit is now r <ti R, which leads to the solution we found 
in previous sections. 

The solution (45) has an immediate interpretation as a D3 brane in a pp wave back- 
ground. In fact, in [38] a class of solutions were found that includes (45) in the case with 
no B2 . Similar solutions were discussed in [39] , although these had a wave profile along the 
D3 brane directions. In the case we are studying, the pp wave is travelling along the D3 
brane, but has a profile orthogonal to the D3 brane. The connection to D3 branes in a pp 
wave background can be made especially explicit in the case of X5 = S^, as the geometry 
can be rewritten as 

d^HB = 777^ {-2dx+dx- + dx" - [aHij + aij\yy{dx+f) + H{rf/^dy^ . (51) 
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Here y are the six coordinates transverse to the D3 brane worldvolume, = t/^ , and aij is 

a tracefree matrix of real numbers. We used the fact that the 1 = 2 harmonics on are 
inherited from quadratic polynomials in M^. This spacetime asymptotes to a homogeneous 
ten dimensional plane wave. The near horizon Schrodingcr group can be thought of as 
a remnant of the large symmetry of the plane wave. The tracefree part of the waveform 
originates from a purely metric ten dimensional solution, whereas the trace is supported by 
the null three form flux. 

6.1 Supersymmetries 

Given that the -D3 brane and the homogeneous pp wave (without flux, B2 = 0) both preserve 
half of the supersymmetries of IIB supergravity, it is natural to expect that this intersection 
of the two solutions will preserve one quarter. 

It is straightforward to check that our solution (9) with S2 = indeed preserves a 
quarter of the supersymmetries. In the absence of NSNS three form flux, all that is needed 
is that the gravitino variation vanish 

Sil^M = {dm + Y^r^-f'^^^rMi^JVPQiis) e = . (52) 
We can introduce the tangent space gamma matrices through 

7" = e"Mr^^ , (53) 
where e^Af is the vielbein. The spinorial covariant derivative is then 

Dm = dM + 7a/3 • (54) 

We now proceed to separate the spinor e into various components using projection 
operators. Firstly one can impose the chirality condition 

7iie = -e , (55) 
where as usual 711 = —7"' i^r- - ^hen we split 

e = e+ + e_ = ^(1 + i7+"^^)e + ^(1 - i7+"^^)e , (56) 

and furthermore split 

e± = el+^ + ei") = -i7"7+e± - h^7"e± ■ (57) 
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Using this decomposition it is straightforward to solve the KiUing spinor equations (52) for 
the near horizon Schrodinger geometry (9). The result is 

4+)=0, e^^=0, e^~^ = V^h{Xr,)n^-\ (58) 

where r]^ ^ is a constant spinor consistent with the projections above and h{X^) is an 
operator that projects onto Killing spinors on the internal space X5. 

Introducing the wave therefore breaks any pre-existing symmetries down to one quarter. 
In particular, if X5 = the background (9) is 1/4 BPS, while if X5 = T^'^ or any other 
Sasaki-Einstein space, the background is 1/16 BPS. The surviving symmetries are half of the 
supertranslations. This suggests that the supersymmetry algebra of this background may 
be the algebra with eight supercharges found in [40]. Other recent interesting discussions 
of a super Schrodinger symmetry for AdS/CFT include [41, 42]. A systematic approach to 
Schrodinger superalgebras was presented in [43]. 

Including a nonzero B2 is likely to break the supersymmetry altogether, as noted in 
[19, 20]. 



7 A renormalisation group flow 

In this section we study a renormalisation group flow. An interesting novelty compared 
with the well studied Poincare invariant renormalisation group flows in AdS/CFT is that 
here the time and spatial parts of the metric can scale differently as a function of the radial 
distance (this fact has been used in [44] and [45] to find flows that interpolate between fixed 
points with different values of the speed of light and z, respectively). This allows more 
freedom in finding solutions. 

The following relatively simple ansatz is consistent: 

dsjiB = { -2s{r)dx+dx- + t{r)dx'^ + ^ - f{X^)m{r){dx+f + p{r)ds\^ \ ,(59) 
F5 = 4i?4(l + ^)vol(X5). (60) 

As before — V^^/ = 12/. The Einstein equations become a set of nonlinear equations for 
the four functions m,p}. There is a first order constraint 

^Jl ±-n 
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and three second order equations 



// I's' t' 5p' 1\ , / s'2 12 8 , 



s t 2 p r J \ r^p r^p^ 

, t' hp' l\ 8s 



t" = {-T-o'z-z)i (64) 



s' bp' 1\ , 8t 

--o r + —5 

s 2 p r J r'^p^ 

Prom these equations and the first order equation (61) it is straightforward to obtain an 
equation for p" . Note the structure of the equations: there are three equations for {s,t,p} 
and then a decoupled equation for m. Therefore any solution will be a (generally non- 
Poincare invariant) renormalisation group flow of the A/" = 4 theory that is then deformed 
into a nonrelativistic background, without backreaction on the initial metric functions. 

Remarkably, it is possible to find an explicit analytic solution to these equations. In 
fact, one can find the general solution with a constant size S^, that is p = 1. The remaining 
radial functions become 



j,4 


4 




+ 4 




4 




-r* 




+ 4 




+ 4 


j,4 





m{r) = [■-^A ^^-^^(1-5), (65) 



8 _ ^8 



r° — r. 



sir) = (l^) ^^7^' (66) 
tir) = (V-:^) ^^V^- (67) 



Here E{k) is a complete elliptic integral of the second kind^. Note that E{0) = 7r/2. There 
is a second independent solution which has a logarithmic divergence as r — > ro, which we 
have not written here. 

The first fact we can note is that this metric asymptotes as r — oo to our original 
solution (9). Thus, we asymptotically develop a Schrodinger symmetry. On the other hand, 
the solution is singular at r = ro- It is easy to check that the Riemann tensor squared 
diverges as r — tq. However, it is a null singularity (this is seen from the fact that it 
takes an infinite boundary time for an infalling null observer to reach the singularity) and 
therefore likely to be admissible as a solution in the full string theory [46]. 

If we had taken m = 0, this flow would have a simple interpretation. The functions 
s and t appearing in the metric (59) are dual to components of the energy momentum 



^Recall that E{k) = f^J^^ \/l — fc^ sin'^ 6d6. The elhptic integral may also be expressed in terms of 
hypcrgcomctric functions or associated Logondro functions. There seems to be a bug in (some versions at 
least of) Maple in simplifying these elliptic integrals. 
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tensor in the (3+1 dimensional) J\f = 4 theory. There is a vacuum expectation value but 

no source for these modes. Therefore we would be looking at the theory in a state in 
which the expectation value for the pressure is anisotropic. This is not especially natural 
in the absence of anisotropic sources. However, once we have m 7^ there is a different 
interpretation. 

Viewed as a dual to a nonrelativistic theory, we identify the x~ coordinate which is now 
no longer a spacetime coordinate in the dual (2+1 dimensional) field theory. According 
to [14, 21], the coefficient of (dx^)"^ couples to the particle number density of the dual 
theory while in our solution the subleading term in m(r) is a constant. Therefore it seems 
that this solution corresponds to a finite energy density but zero particle density. This 
statement is furthermore supported by the absence of a {dx~)'^ term, which appeared in 
the finite particle density solutions of [19, 20, 21]. 

Another consistent simplification of the equations above is to set s = t. We have not 
been able to solve the resulting equations analytically. This flow, with m = 0, would be a 
Poincare invariant state of the M = 4 theory and therefore less specific to the nonrelativistic 
setup. 

8 Coulomb branch solutions 

The methods we have been using are easily adapted to produce Coulomb branch solutions 
that asymptote to nonrelativistic backgrounds. Our starting point will be the explicit 
supersymmetric and asymptotically AdS backgrounds that were constructed in [47] . There 
are five such backgrounds, which preserve an SO{n) x SO(6 — n) subgroup of the SO{&j 
symmetry of the M = 4 theory. The fact that these backgrounds only involve the metric and 
five form field strength make them especially amenable to non-relativistic generalisation. 
Some solutions related to those we shall present shortly, but with different asymptotics, 
were found in [48]. 
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8.1 Asymptotically Schrodinger solution with S0(2) x 50(4) symmetry 

We will write down directly the Coulomb branch solution with the nonrelativistic deformation.^ 
The metric and five form of our solution are 

„2 



ds^ = Cr^ (^-2dx+dx- + dx"^ + - m(r, S^){dx+)^^ 



^ ,^2 cos^ 9 sin^ 9 ,^2 
+ Cd9^ + —^dn^a + dn^i , 

F5 = [r (L^ + 4r^ + cos 29) dr - L^r^ sin 29d9] A dx+ A dx~ AdxAdy, (68) 



with the remaining components of F5 obtained by Hodge dualising. We introduced the 
functions 

^^ = 1 + 7^' C^ = l + 7^cos20. (69) 

Here L gives a notion of the distance between the smeared D'i branes. At L = the 
solution reduces to a gravitational wave in pure AdS. Setting m = gives the solution in 
[47]. Introducing a nonzero m will still solve the HB equations of motion (11) provided that 
m is harmonic on the spacetime with m = 0, i.e. V^m = 0. Writing this out explicitly gives 

(70) 

The differential equation (70) can be separated easily if we take m to be independent of 
the coordinates of the and S"^ in the spacetime (68). To find a regular solution on the 

we must therefore find spherical harmonics on that only depend on the coordinate 9. 
This is the coordinate that arises if we parametrise C by 

X = cos 9n4 + sin 9 n2 , (71) 

where n4 and n2 are unit vectors taking values in an orthogonal and M^, respectively, in 
M^. Recall that spherical harmonics on with eigenvalue £{£ + 4) are given by tracefree 
symmetric tensors of rank £ on R^. Requiring 50(4) x 50(2) invariance restricts us to 
tensors of the form 

M 

c„,..„,^x«i • • • = J2 «Ar(x? + xl + xl + xD^'ixl + x^)^-^ . (72) 

N=0 

The tracefree condition (after allowing for the symmetry in this tensor) then imposes M 
relations between the M + 1 different a^- Thus we can find a spherical harmonic solution 



^In this section we set the AdS radius J? = 1 for simplicity. 



18 



for any even I = 2M. By our formula (40) we thus obtain Coulomb branch solutions that 

asymptote to nonrelativistic backgrounds for any positive integer z € Z+. Let us write 
explicitly the solution in the case i = 2. Using the parametrisation (71) and furthermore 
solving the remaining radial equation in (70) leads to the solution 

m(,=2) =a(^l + 1^ j (^-^+ cos^ . (73) 

Here a is a free parameter. There may of course be other more complicated solutions to 
(70) in which the variables do not separate. It is pleasing that the simplest regular solution 
turns out to have the properties we are interested in. The solution is a gravitational wave 
in a supersymmetric Coulomb branch geometry that asymptotes to our nonrelativistic dual 
spacetime (9), with a particular choice for the £ = 2 spherical harmonic in that solution. 

8.2 S0{3) X SO{3) and SO{5) symmetric solutions 

The treatment of the other two possible symmetries is entirely analogous, and so we shall 
be brief. We take the SO (3) x SO (3) case first. The solution is 

ds"^ = Cr^Xl-2dx+dx- + dx'^ + ^-m{r,S^)idx+)^\ 



+ jde^ + -^^dnl2 + — ^ — dnl^ , (74) 

+ SLV^ + Sr^ + L2(l2 + 2r2) cos 29 



F5 = 



-dr 



2VL2 + r2 

- L\\/L^~+r^ sin 29 dO A dx+ A dx~ A dx A dy , (75) 



with the remaining components of F5 obtained by Hodge dualising. While C is as before in 
(69), we now have 

X' = l + ^. (76) 

As before we can find solutions with integer values of z. The solution that has an asymptotic 
Schrodinger symmetry is 



'^^'='^ = " LWL^ + r^ - + cos ^ ) , (77) 



where again a is a free parameter. 
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For the case with 50(5) symmetry, the solutions take the form 



ds'' = ^ (^-2dx+dx- + dx^ + ^- m{r, S^){dx+f^ + (A^d^^ + cos^ ^^^^^ ^ 



r\3L^ + l2L'r' + Sr^ + L^(3L^ + 2r') cos 20) 



2(L2 + r2)3/2 



dr 



2Ji 



:sin26'd6' 



A dx'^ Adx Adx Ady , 



(78) 



with the remaining components of obtained by Hodge duahsing. For this solution we 
have 



A 



12 



1 + 



L2 



The solution with an asymptotic Schrodinger symmetry has 



-- +COS 

6 



(79) 



(80) 



Unlike in the other cases, there are two parameters here, a and b. Both of these solutions 
have an asymptotic Schrodinger symmetry. One of them has a divergence as r ^ 0, but 
given that the background itself is singular in this limit it is not clear that this is sufficient 
a pathology to throw away the solution. As before, there are also solutions for all other 
positive integer z. 

There were two further metrics considered in [47]. These are related to the 50(2) x 50(4) 
and 50(5) metrics presented above by letting —L?'. The metric with 50(5) symmetry 

and the sign of inverted is believed to be nonphysical. 

We found above that introducing our pp wave into AdSr, x broke all the supersym- 
metries except for half of the supertranslations. The Coulomb branch solutions for A/' = 4 
super Yang Mills theory break the superconformal symmetries but preserve all the super- 
translations [47]. Therefore, it seems likely that the solutions we have found in this section 
will preserve half of the supertranslation symmetries. 



9 Discussion 

In this paper we have presented several solutions to IIB supergravity with features that 
seem likely to be of interest as duals for nonrelativistic quantum critical theories in 2+1 
dimensions. We initiated a study of basic features of these solutions including stability, 
regularity and renormalisation flow generalisations. In this discussion we note important 
open questions and comment on the dual interpretation of these backgrounds. 
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Given that our backgrounds have (at least) a manifest 2+1 dimensional Galilean symme- 
try, it seems clear that they are dual to nonrelativistic theories defined on a 2+1 dimensional 
spacetime. In particular, the mixing between the 'AdS' and internal coordinates in the bulk 
should not alter this fact. Furthermore, our z = 2 solutions arc continuously connected to 
the embedding found in [19, 20, 21]. It therefore seems likely that, as with those solutions, 
our bakgrounds are dual to (differing) twisted DLCQs of = 4 Super Yang-Mills theory. 
The dual description of the solutions with 2; 7^ 2 is possibly less clear from this point of 
view. However, we noted that they share with the z = 2 cases the interpretation as the 
near horizon limit of a branewave. 

It is of interest for several reasons to generalise the solutions we have found to a finite 
temperature and also a finite number density, with a (dx~)^ term in the metric. One rea- 
son is that most potential applications to experimental systems will be at least at finite 
number density and often at finite temperature. Another reason is that, as we have men- 
tioned several times above, in solutions with vanishing temperature and number density 
the periodically identified circle (3) is null. This is expected to invalidate the supergravity 
approximation as wrapped strings can become light [20]. 

The mixing in our solution (9) between the noncompact and the X5 directions makes it 
hard to find generalisations with a (dx")^ component in the metric. This is what we would 
need to make the identified circle spacelike. In fact, all of the solutions discussed in this 
paper seem to be an application of the Garfinkle-Vachaspati solution generating technique 
[49, 50] (for a recent review of this method in a string theory context see for instance [38]). 
This method requires a null (hyper surface-orthogonal) Killing vector field, which will not 
be present in finite temperature solutions, or more generally if d^-- is spacelike. A different 
approach will be required to find the finite temperature solutions. 

Restricting to the solutions we have found in this paper, there are several immediate 
remaining questions. We did not analyze the stability of the solutions with z ^ 2 and 
we did not consider the stability of any of the backgrounds against other modes such as 
gravitational perturbations. We also did not find the explicit form of the unstable scalar 
mode. 

It seems clear that our constructions can straightforwardly be generalised to different 
spacetime dimensions. It would also be interesting to generalise the renormalisation group 
flow solution of section 7 to include a B2 field, and also to different values of z. 
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